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Abstract 

In this paper, the performance of a binary phase shift keyed random time-hopping impulse radio system with 
pulse-based polarity randomization is analyzed. Transmission over frequency-selective channels is considered and 
the effects of inter-frame interference and multiple access interference on the performance of a generic Rake receiver 
are investigated for both synchronous and asynchronous systems. Closed form (approximate) expressions for the 
probability of error that are valid for various Rake combining schemes are derived. The asynchronous system 
is modelled as a chip-synchronous system with uniformly distributed timing jitter for the transmitted pulses of 
interfering users. This model allows the analytical technique developed for the synchronous case to be extended to 
the asynchronous case. An approximate closed-form expression for the probability of bit error, expressed in terms of 
the autocorrelation function of the transmitted pulse, is derived for the asynchronous case. Then, transmission over 
an additive white Gaussian noise channel is studied as a special case, and the effects of multiple-access interference 
is investigated for both synchronous and asynchronous systems. The analysis shows that the chip-synchronous 
assumption can result in over-estimating the error probability, and the degree of over-estimation mainly depends on 
the autocorrelation function of the ultra-wideband pulse and the signal-to-interference-plus-noise-ratio of the system. 
Simulations studies support the approximate analysis. 
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I. Introduction 

Since the US Federal Communications Commission (FCC) approved the limited use of ultra-wideband (UWB) 
technology [1], communications systems that employ UWB signals have drawn considerable attention. A UWB 
signal is defined to possess an absolute bandwidth larger than 500MHz or a relative bandwidth larger than 20% and 
can coexist with incumbent systems in the same frequency range due to its large spreading factor and low power 
spectral density. UWB technology holds great promise for a variety of applications such as short-range high-speed 
data transmission and precise location estimation. 

Commonly, impulse radio (IR) systems, which transmit very short pulses with a low duty cycle, are employed to 
implement UWB systems ([2]-[6]). In an IR system, a train of pulses is sent and information is usually conveyed 
by the position or the polarity of the pulses, which correspond to Pulse Position Modulation (PPM) and Binary 
Phase Shift Keying (BPSK) 5 , respectively. In order to prevent catastrophic collisions among different users and 
thus provide robustness against multiple-access interference, each information symbol is represented by a sequence 
of pulses; the positions of the pulses within that sequence are determined by a pseudo-random time-hopping (TH) 
sequence specific to each user [2]. The number Nf of pulses representing one information symbol can also be 
interpreted as pulse combining gain. 

In "classical" impulse radio, the polarity of those Nf pulses representing an information symbol is always the 
same, whether PPM or BPSK is employed ([2], [7]). Recently, pulse-based polarity randomization was proposed, 
where each pulse has a random polarity code (±1) in addition to the modulation scheme ([8], [9]). The use of 
polarity codes can provide additional robustness against multiple-access interference [8] and help optimize the 
spectral shape according to FCC specifications by eliminating the spectral lines that are inherent in IR systems 
without polarity randomization [10]. 

A TH-IR system with pulse-based polarity randomization can be considered as a random CDMA (RCDMA) 
system with a generalized signature sequence, where the elements of the sequence take values from {— 1,0, +1} 
and are not necessarily independent and identically distributed (i.i.d.) [8]. The performance of RCDMA systems with 
i.i.d. binary spreading codes has been investigated thoroughly in the past (see e.g. [11]-[13]). Recently, [14] and [15] 

have considered the problem of designing ternary codes for TH-IR systems. Moreover, in [8], the performance of 

5 Since IR is a carrierless system, the only admissible phases are and it. Therefore, BPSK becomes identical to Binary Amplitude-Shift 
Keying (BASK) in this case. 
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random TH-IR systems with pulse-based polarity randomization has been investigated over additive white Gaussian 
noise (AWGN) channels, assuming symbol-synchronized users. To the best of our knowledge, no study concerning 
the bit error probability (BEP) performance of Rake receivers (with various combining schemes) for a random TH-IR 
system with pulse-based polarity randomization in a multiuser, frequency-selective environment has been reported 
in the literature. In this paper, we investigate such a system and provide (approximate) closed-form expressions for 
its performance. We consider an important case in practice, where the different users are completely asynchronous. 
We begin by considering the chip-synchronous case where the symbols of different users are misaligned but this 
misalignment is an integer multiple of the chip interval. Subsequently, we treat a more general asynchronous 
case, where we show that the system can be represented as a chip-synchronous system with uniform timing jitter 
between zero and the chip interval for each interfering user. We consider frequency-selective channels and analyze 
the performance of Rake receivers with various combining schemes. 

The remainder of the paper is organized as follows. Section II describes the transmitted signal model for a 
TH-IR system with pulse-based polarity randomization. In Section III, both chip-synchronous and asynchronous 
systems over frequency-selective channels are considered, and the performance of Rake receivers is analyzed for 
various combining schemes. Simulation studies are presented in Section IV, followed by some concluding remarks 
in Section V. 

II. Signal Model 

We consider a BPSK random TH-IR system with N u users, where the transmitted signal from user k is represented 

by 



where wt x {t) is the transmitted UWB pulse with duration T c , Ef. is the bit energy of user k, Tf is the "frame" 
time, Nf is the number of pulses representing one information symbol, and frj^/^j e {+1, — 1} is the information 
symbol transmitted by user k. In order to allow the channel to be shared by many users without causing catastrophic 



are independent for (k, j) / This TH sequence provides an additional time shift of c- T c seconds to the jth 




(1) 



collisions, a time-hopping sequence {c- } is assigned to each user, where c- G {0, 1,...,N C — 1} with equal 
probability, with N c denoting the number of possible pulse positions in a frame (N c = Tf/T c ), and and cf^ 
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Fig. 1. A TH-IR signal with pulse-based polarity randomization where Nf — 6, N c — 4 and the TH sequence is {2, 1, 2, 3, 1, 0}. Assuming 
that +1 is currently being transmitted, the polarity codes for the pulses are {+1, +1, — 1, +1, — 1, +1}, 



pulse of the kth user where the pulse width T c is also considered as the chip interval. 

N = NfN c represents the total processing gain of the system. Due to the regulations by the FCC [1], each user 
can transmit a certain amount of energy in a given time interval. Since the symbol (bit) energy of the signal defined 
in (0 is constant (denoted by E^), we consider a fixed symbol interval; hence, a constant total processing gain N 
throughout the paper. 

(k) 

The random polarity codes d - s are binary random variables taking ±1 with equal probability, and such that 
dj and of are independent for (k,j) ^ [8]. Use of random polarity codes helps reduce the spectral lines 
in the power spectral density of the transmitted signal [10] and mitigate the effects of MAI [8]. The receiver for 
user k is assumed to know its polarity code. 

Defining a sequence {Sj} as 



(k) 



Otherwise. 



we can express (0 as 




i fc) w = \ht E *V ^V/N^ut - jT c ), o) 



J j=—oo 



which indicates that a TH-IR system with polarity randomization can be regarded as an RCDMA system with a 
generalized spreading sequence {sj} ([16], [8]). Note that the main difference of the signal model in (0 from the 
"classical" RCDMA model ([11]-[13]) is the use of {— 1,0, +1} as the spreading sequence, instead of {— 1,+1}. 
The system model given by equation (0 can represent an RCDMA system with a processing gain of Nf, by 
considering the special case when Tf = T c . 

An example TH-IR signal is shown in Figure where six pulses are transmitted for each information symbol 



(Nf = 6) with the TH sequence {2, 1, 2, 3, 1, 0}. 

III. Performance Analysis 
We consider transmission over frequency selective channels, where the channel for user k is modelled as 

L 

h^(t) = ^ k) 5(t -(I- 1)T C - r k ), (4) 
l=i 

(k) 

where a, and are the fading coefficient of the Ith path and the delay of user k, respectively; T c is the minimum 
resolvable path interval. We set t\ = without loss of generality. We assume that the channel characteristics remain 
unchanged over a number of symbol intervals, which can be justified by considering that the symbol duration in a 
typical application is on the order of tens or hundreds of nanoseconds, and the coherence time of an indoor wireless 
channel is on the order of tens of milliseconds. 

Note that the channel model in © is quite general in that it can model any channel of the form Ylb=i o[ k ^5(t— r, ) 
if the channel is bandlimited to 1/T C . Thus, each realization of an arbitrary (and nonuniformly sampled) channel 
model, e.g., the 802.15.3a UWB channel model [17], can be represented in the form of equation (4). Only the 
statistics of the tap amplitude are changed when the tap spacing is changed to a uniform spacing. 

Using © and ©, the received signal can be expressed as follows: 

n u rrr oo 

r ^ = E\ W £ dfb^ NA 4 k \t-jT f -cfT c -r h ) + a n n{tl (5) 

k=l V ' j=—oo 

where n(t) is a white Gaussian noise with zero mean and unit spectral density, and 

L 

) (t) = ^2a { l k) w rx (t-(l-l)T c ), (6) 



it 

l=i 

with w rx (t) being the received UWB pulse with unit energy. 

We consider a Rake receiver for the user of interest, say user 1, and express the template signal at the Rake 
receiver as follows: 



(1+1)^-1 

»2Lp(*)= E dfv{t-jT f -cfT c \ (7) 
j=iN f 



where 

L 

V(t) =Y,PlWrx(t-{l-l)T c ), (8) 
1=1 

with /3 = \Px, Pl] being the Rake combining weights. 

The template signal given by (Q and (JSJl can represent different multipath diversity combining schemes by 
choosing an appropriate weighting vector (3: In an M-finger Rake the weights for (L — M) multipath components 
not used in the Rake receiver are set to zero while the remaining M weights are determined according to the 
combining scheme, such as "Equal Gain Combining (EGC)" or "Maximum Ratio Combining (MRC)". 

The output of the Rake receiver can be obtained from ©-(JSjl as follows: 

L 

yi = bf ] y/Ejff Y,a\ 1] Pi + a + a + n, (9) 
1=1 

where the first term is due to the desired signal, a is the self interference of the received signal from user 1 itself, 
which we call inter-frame interference (IFI), a is the MAI from other users and n is the output noise, which is 
approximately distributed as n ~ N (o , N f al YjL x 0f) for large N f (Appendix EJ. 

Inter-frame interference (IFI) occurs when a pulse of user 1 in a frame spills over to an adjacent frame due 
to the multipath effect and consequently interferes with the pulse in that frame (Figure 12). The IFI in © can be 
expressed, from (J5jl and Q, as 

r^-(i+l)N f -l 

& = \W f ^ a ™' (10) 

where 

oo 



a m = d« Y, 4 %/N fi ^ (C? " m)T, + {cf - c«)T c ) , (11) 



j=-oo 



with 4>uv{x) denoting the cross-correlation between u^ k \t) of © and v(t) of ©: 

POO 

^){x)= / u^{t-x)v{t)dt. (12) 



Note that a m in (fTXt denotes the IFI due to the transmitted pulse in the mth frame of user 1, and the sum of 
such IFI terms over Nf frames is equal to a, as seen in (fTQl) . In Appendix |Bj we show that these Nf terms form a 
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Fig. 2. Inter-frame interference (IFI) from the (m — l)th frame to the mth frame, where p„ denotes the position of the first user's pulse 
in the mth frame. Only the signals from the mth frame of the template (the signal on the top) and from the (m — l)th frame of the first 
user are shown. The IFI can also result from a spill-over of the signal at the mth frame of the template to the (m + l)th frame when any 
of the pulses of the first user in the (m + l)th frame overlap with those pulses that spill over. 



1 -dependent sequence 6 when L < N c + 1 and their sum converges to a Gaussian random variable for a large Nf. 
This result is summarized in the following lemma: 

Lemma 3.1: As Nf — > oo, the IFI a in ( fTUl is asymptotically normally distributed as 



n L - 1 



a ~ 



c j=i li=i 



L-j 



2\ 



(13) 



for L < N c + 1. 

Proof: See Appendix |B] 

Note that for a Rake receiver with one finger such that 0\ = 1 and /3j = for / = 2, L, the expression reduces 

toa~jv(o,^££i l i(~ (1) ^ 2 



AT C 2 Z^/=l 1 V a i+i 

Due to the FCC's regulation on peak to average ratio (PAR), Nf cannot be chosen very small in practice. Since 
we transmit a certain amount of energy in a constant symbol interval, as Nf gets smaller, the signal becomes more 
peaky as shown in Figure |3] Therefore, the approximation for large Nf values can be quite accurate for real systems 
depending on the system parameters. 

When L > N c + 1, the pulses in a frame always spill over to the adjacent frame(s). In this case, the Nf terms in 
(fTOl form a \(L— l)/N c ] -dependent sequence and the asymptotic distribution of the IFI is given by the following 
lemma: 

6 A sequence {X n } ne z is called a D-dependent sequence, if all finite dimensional marginals (X ni , X ni ) and (X mi , X m , ) are 
independent whenever mi — m > D. 
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Fig. 3. Two different cases for a BPSK-modulated TH-IR system with pulse-based polarity randomization when N = 24. For the first 
case, N c = 8, Nf — 3 and the pulse energy is £7/3; for the second case, N c — 4, iV/ = 6 and the pulse energy is £7/6. 



Lemma 3.2: As Nf — > oo, the IFI a in dlOI > is asymptotically normally distributed as 

rx-j 



a ~ 



L-7V C 



J'=l 



,z=i 



r - -i 



(14) 



for L > N c + 1. 

Proof: See Appendix ICl 

The MAI term in © can be considered as the sum of MAI terms from each user, that is, a = Ylk=2 where 
each a( fc ) is in turn the sum of interference due to the signals in the frames of the template: 

■(i+l)JV/-l 



*(*) 




J m=iNf 



with 



(15) 



® = <® E 4X/^j^ ) (^- m ) r /+(4* ) -^ ) ) r «+ 7 *)' 

j=-oo 



(16) 



where cjJuv (x) is as in (II 2b and is the delay of the fcth user. 

The effects of MAI will be different for synchronous and asynchronous systems, as investigated in the following 
subsections. 



A. Symbol-Synchronous and Chip-Synchronous Cases 

In the symbol-synchronous case, the symbols from different users are exactly aligned. In other words, = 
for k = 2, ...,N U . On the other hand, for a chip-synchronous scenario, the symbols are misaligned but the amount 



of misalignment is an integer multiple of the chip interval T c . That is, r& = AfcT c , for k = 2, iV M , where A& is 
uniformly distributed in {0, 1, N — 1} with N = N c N f . 

Note that the assumption of synchronism may not be very realistic for a UWB system due to its high time 
resolution. However, the aim of this subsection is two-fold. First, we will show that the BEP performance of 
the UWB system is the same whether the users are symbol-synchronized or chip-synchronized. Second, we will 
extend the result for the chip-synchronous case to a more practical asynchronous case by modelling asynchronous 
interfering users as chip-synchronous users with uniform timing jitter, as will be shown in the next subsection. 

The following lemma gives the asymptotic distribution of MAI from a user for a large number of pulses per 
symbol. 

Lemma 3.3: As Nj — > oo, the MAI from user k, which is chip-synchronized to user 1, is asymptotically 
normally distributed as 



L f j \ 2 L-1 / j - 21 



j=X \i=i J j=i \i=i 



(17) 



The result is also valid for a symbol-synchronous scenario. 
Proof: See Appendix iDl 

Note that when (3\ = 1 and 0i = 0, for / = 2, L, we have ~ N ^0 , jf- Yla=i a / 2 )> which represents the 
result for a Rake receiver with a single finger that picks up the first path signal component only. 

Note that the Gaussian approximation in Lemma 1331 is different from the standard Gaussian approximation (SGA) 
used in analyzing a system with many users ([19]-[21]). Lemma 1331 states that when the number of pulses per 
information symbol is large, the MAI from an interfering user is approximately distributed as a Gaussian random 
variable. 

We also note from Lemma 13.31 that the effect of the MAI is the same for symbol-synchronized and chip- 
synchronized cases. This is due mainly to the pulse-based polarity randomization, which makes the probability 
distribution of the MAI independent of the information bits of the interfering user, as can be observed from d 1 6b . 
Since the probability that a pulse of the template signal overlaps with any of the pulses of an interfering user is 
the same whether the users are symbol-synchronous or chip-synchronous, the probability distributions turn out to 
be the same for both cases. 

An approximate expression for BEP can be derived from using Lemma l3~Tl Lemma l3~2l and Lemma l3~3l as 
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follows: 



Pe-Q 



N C N ( 



a IFI,l + JV a IFI,2 + iV Efc=2 Ek&MAI,!: + 2-4=1 A 2 



(18) 



where 



2 

a IFIA 



2 

a IFI,2 



L-j 



j=i 



min{iV c ,L}-l 
3=1 



0. 



(19) 



(20) 



and 



= E f E A^il + E f E «i fc) a+w 

j=l \i=l / j=l \i=l / 



(21) 



Equation (fTHl implies that, for a fixed total processing gain N, increasing N c , the number of chips per frame, 
will decrease the effects of IFI, while the dependency of the expressions on the MAI remains unchanged. Hence, 
an RCDMA system, where Nf = N, can suffer from IFI more than any other TH-IR system with pulse-based 
polarity randomization, where Nf < N, if the amount of IFI is comparable to the MAI and thermal noise. 

B. Asynchronous Case 

Now consider a completely asynchronous scenario. In this case, it is assumed that in (II 6b is uniformly 
distributed according to U[0, NT C ) for k = 2, N u . 

In order to calculate the statistics of the MAI term in ®, the following simple result will be used. 

Proposition 3.1: The MAI in the asynchronous case has the same distribution as the MAI in the chip-synchronous 
case with interfering user k having a jitter e^, for k = 2, N u , which is the same for all pulses of that user and 
is drawn from the uniform distribution U[0,T c ). 

Proof: Consider dl6b . For k = 2, . . . , N u , is uniformly distributed in the discrete set {0, T c , (N— 1)T C } in the 
chip-synchronous case. In the asynchronous case, t^. is a continuous random variable with distribution li[0, NT C ). If 
the jitter in the chip-synchronous case is uniformly distributed with U[0, T c ), then T^ + e^ is uniformly distributed 
as li[0, NT C ) hence is equivalent to the distribution of t^. in the asynchronous case. □ 

Proposition 13.11 reduces the performance analysis of asynchronous systems to the calculation of the statistical 
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properties of 

oo 

«L fc) = d® £ d f )h WN^ (U - m)T, + (cf - c«)T c + r k + e fc ) , (22) 

J=-oo 

where = A k T c takes on the values {0, T c , (N — 1)T C } with equal probabilities and e k ~ W[0, T c ). This 
problem is similar to the analysis of TH-IR systems in the presence of timing jitter, which is studied in [18]. 
However, in the present case, the timing jitter of all pulses of an interfering user is the same instead of being i.i.d. 

The following lemma approximates the distribution of the MAI from an asynchronous user, conditioned on the 
timing jitter of that user when the number of pulses per symbol, Nf, is large. 

Lemma 3.4: As Nf — ► oo, the MAI from user k given e k has the following asymptotic distribution: 



a< fc > 



|e*~Af(0, ^* 2 MAI , k (e k )\ , (23) 



where 

L-l / j 



2 



°MAi,k{tk) = £ EAK+Li-i-R^ " + + (3 3+ iafR{T c - e k ) 

j=0 \i=i ) 

L-l / j \ 2 

+ £ ^ \Pl+L-j-iR(e k ) + (3 l+L ^R(T c - e k )] + afl x (3 L R(e k ) , (24) 
j=0 \i=i ) 

with R(x) = f™^ w TX {t - x)w rx {t)dt. 
Proof: See Appendix IE] 

Note that when e k = 0, which corresponds to the chip-synchronized case, (l2"4l reduces to (l2"Tl i. 

From Lemma l3~4l we can calculate, for large Nf, an approximate conditional BEP given e = [e2 ■ ■ ■ £jyj as 

(1) 

~XJV a IFI,l + W a IFI,2 + ]?" Sfc=2 EkO~ 2 MAI k {e k ) + <T 2 X)z=l A" 



where o\ lAI k (e k ) is as in d2"4l and <r| F/ x and <t| fj 2 are as in dT^l and respectively. 

By taking the expectation of (l2"5t with respect to e = [e2, . . . , ejv„]> where ~ W[0, T c ) for /c = 2, N u , we 
find the BEP: 

1 /" Tc /" T ° , 

TV -1 / - • • / 1 6 ^ e2 ' - - de N u - 

T c " Jo jo 
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However, when the number of users is large, calculation of (|2"6T ) becomes cumbersome since it requires integration 
of P e \e over (N u — 1) variables. In this case, the SGA [19]-[21] can be employed in order to approximate the BEP 
in the case of large number of equal energy interferers: 

Lemma 3.5: Assume that all the interfering users have the same bit energy E. Then, as N u — > oo, a/y/N u — 1, 
where a is the MAI term in ©, is asymptotically normally distributed as 

where 

L-l Tc / j \ 2 

H°MAI,k{tk)} = 7f E / E WS^-i^c - e fc ) + af^Rieu)] + ^ j+1 af R{T C - e h ) de k 
Ic j= o Jo \i=i J 

l L-l Tc / j \ 2 

+ 7^E / E a ! fc) [A+L-i-i^lefc) + Pi+L-jRiTc - e k )} + af + ) 1 /3 Lj R(e fe ) de fc . (28) 
lc j= o Jo \i=i J 
Proof: See Appendix |F| 

The BEP can be approximated from Lemma 1331 as 

p^ Q l ^_ ^ Sfel S ] (29) 

V V 5v5^W + 5 a ^,2 + §(N U - l)E{al IAI k (e k )} + al Ef =1 Pf 

for large A^j and N u , and for equal energy interferers. 



From ( I29l l we make the same observations as in the synchronous case. Namely, for a given value of the total 
processing gain N = N c Nf, the effect of the MAI on the BEP remains unchanged while the effect of the IFI 
increases as the number of chips per frame, N c , decreases. Hence, the IFI could be more effective for an RCDMA 
system, where N c = 1. 

C. Different Rake Receiver Structures 

In the previous derivations, we have considered a Rake receiver with L fingers, one at each resolvable multipath 
component (see © and (JSJl). A Rake receiver combining all the paths of the incoming signal is called an all-Rake 
(ARake) receiver. Since a UWB signal has a very large bandwidth, the number of resolvable multipath components 
is usually very large. Hence, an ARake receiver is not implemented in practice due to its complexity. However, it 
serves as a benchmark for the performance of more practical Rake receivers. A feasible implementation of diversity 
combining can be obtained by a selective-Rake (SRake) receiver, which combines the M best, out of L, multipath 



13 

components. Although an SRake receiver is less complex than an ARake receiver, it needs to keep track of all 
the multipath components and choose the best subset of multipath components before feeding it to the combining 
stage. A simpler Rake receiver, which combines the first M paths of the incoming signal, is called a partial-Rake 
(PRake) receiver [22]. 

The BEP expressions derived in the previous subsections for synchronous and asynchronous cases are general 
since one can express different combining schemes by choosing appropriate combining weight vector, (3. For 
example, if we consider the maximum ratio combining (MRC) scheme, the weights can be expressed as follows 
for ARake, SRake and PRake receivers: 

1) ARake: In this case, the combining weights are chosen as (3 = where (3 = [Pi ■ ■ ■ /3l] are the Rake 
combining weights in © and = [ap . . . aP] are the fading coefficients of the channel for user 1. 

2) SRake: An SRake receiver combines the best M paths of the received signal. Let S be the set of indices of 
these best fading coefficients with largest amplitudes. Then, the combining weights /3 in ((SJl are chosen as follows: 



A 



af\ leS 

(30) 

0, 1<£S 



3) PRake: A PRake receiver combines the first M paths of the received signal. Therefore, the weights of an 
SRake receiver with MRC scheme are given by the following: 



01 



af\ l = l,...,M 

(31) 

0, l = M + l,...,L 



where M < L. 

D. Special Case: Transmission over AWGN Channels 

From the analysis of frequency-selective channels, we can obtain the expressions for AWGN channels as a special 
case, which might be useful for intuitive explanations. 

Considering the expressions in ©-(|8l, and setting a\ = [3\ = 1 and ai = fti = for I = 2, ...,L, the output of 
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the matched filter (MF) receiver can be expressed as 

Vl = ^EjTfb^ +a + n, (32) 

where the first term is the signal part of the output, a is the multiple-access interference (MAI) due to other users 
and n is the output noise, distributed as n ~ M(0, Nfafy. Note that there is no IFI in this case since a single path 
channel is assumed. 

The MAI is expressed as a = ^^ 2 a' fc ', where the distribution of in the symbol-synchronous and chip- 
synchronous cases can be obtained from Lemma 13.31 as 



l(fcW (°'|r)- (33) 



Then, the BEP can be obtained as follows: 



L LA —A (34) 

VV N2Zh=2 E h + < J 

where N = N c Nf, which is the total processing gain of the system. Note from (134b that the BEP depends on iV c 
and Nf only through their product. Hence, the system performance does not change by changing the number of 
symbols per information symbol Nf and the number of chips per frame N c as long as N c Nf is held constant. This 
is different from the general case of (fT8b . where the IFI is reduced for larger N c . Therefore, for AWGN channels, 
the BEP performance of a TH-IR system with pulse-based polarity randomization is the same as the special case 
of an RCDMA system. 

Considering [8], the BEP for TH-IR systems without pulse-based polarity randomization is given by the following 
expression for the case of a synchronous environment with a large number of equal energy interferers: 

/ 



Q 



(35) 



(N u -1)§ (l + ^)+o2 

where E is the energy of an interferer. 

Comparing d34t and d35L we observe that, for Nf > 1, the MAI affects a TH-IR system without polarity 
randomization more than it affects a TH-IR system with pulse-based polarity randomization and that the gain 
obtained by polarity randomization increases as Nf increases (in an interference-limited scenario). The main reason 
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behind this is that random polarity codes make each interference term to a pulse of the template signal (see d!6l >) 
a random variable with zero mean since it can be plus or minus interference with equal probability. On the other 
hand, without random polarity codes, the interference terms to the pulses of the template signal have the same sign, 
hence add coherently, which increases the effects of the MAI. 

Note that the effects of the MAI reduce if the UWB system without pulse-based polarity randomization is in 
an asynchronous environment. Because, in such a case, the MAI terms from some of the pulses add up among 
themselves while the remaining ones add up among themselves and the polarities of these two groups are independent 
from each other. Hence, the average MAI is smaller than the symbol-synchronous case but it is still larger than or 
equal to the MAI for the UWB system with pulse-based polarity randomization, where the sign of each interference 
term is independent (see [23] for the trade-off between processing gains in TH-IR systems with and without polarity 
randomization). 

For TH-IR systems with polarity randomization, we can approximate, using Lemma 1331 the total MAI in the 
asynchronous case for a large number of equal energy interferers as 



Let 7 = £ J Q Tc R 2 (e)de = ± J^R 2 (e)de. Then, from 03, a ~ M (0 , y(N v - l)E/N c ). Note from O 
that for equal energy interfering users, the MAI in the symbol/chip-synchronous case is distributed as a ~ 
A/*(0, (N u — l)E/N c ). Hence we see that the difference between the powers of the MAI terms depends on the 
autocorrelation function of the UWB pulse. For example, for the autocorrelation function of d39l > below, 7 0.2 
and symbol/chip-synchronization assumption could possibly result in an over-estimate of the BEP depending on 
the signal-to-interference-pulse-noise ratio (SINR) of the system. 

From d32t and d36L the BEP of an asynchronous system can be approximately expressed as follows: 



for large values of N u . Similar to the synchronous case, the performance is independent of the distribution of N 
between iV c and Nf. Therefore, the TH-IR system performs the same as an RCDMA system in this case. 




(36) 




(37) 
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Fig. 4. UWB pulses and autocorrelation functions for T c = 0.5ns. 

E. Average Bit Error Probability 

In order to calculate the average BEP, the previous expressions for probability of bit error need to be averaged 
over all fading coefficients. That is, P avg = E{P e (a^\ . . . , a^)}, which does not lend itself to simple analytical 
solutions. However, this average can be evaluated numerically, or by Monte-Carlo simulations. 



IV. Simulation Results 

In this section, the BEP performance of a TH-IR system with pulse-based polarity randomization is evaluated by 
conducting simulations in MATLAB. The following two types of (unit energy) UWB pulses and their autocorrelation 
functions are employed as the received UWB pulse w rx {t) in the simulations (Figure 0: 



Wl {t)=[l r2 
Ri(At) 

w 2 {t) 



r or 



1 



0.5T r < t < 0.5T r 



R 2 (At) = < 



-Ai/T c -r-l, < At < T c 
At/T c + l, -T c < At < 



(38) 
(39) 
(40) 

(41) 



where E p of wi(t) is the normalization constant, t = T c /2.5 is used in the simulations, and the rectangular pulse 
W2(t) is chosen as an approximate pulse shape in order to compare the performance of the system with different 
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Symbol/Chip-sync w/ polarity rand, Theo 
Async w^t) w/ polarity rand, Theo 
, Async w 2 (t) w/ polarity rand, Theo 
—t— Symbol-sync w/o polarity rand, Theo 
-©- Symbol-sync w/ polarity rand, Simu 
Chip-Sync w/ polarity rand, Simu 
D Async w^t) w/ polarity rand, Simu 
Async w 2 (t) w/ polarity rand, Simu 
-A- Symbol-sync w/o polarity rand, Simu 



SINR (dB) 

Fig. 5. BEP vs SINR for different cases, where N c = 5, Nf — 15, N u — 10, i?i = 0.5 and E = 1. Transmission over an AWGN channel 
is considered. 



pulse shapes. 

Figure |5] shows the BEP performance of a 10-user system (N u = 10) over an AWGN channel, where Nf = 15 
and N c = 5. The bit energy of the user of interest, user 1, is E% = 0.5, whereas the interfering users transmit bits 
with unit energy (E^ = 1 for k = 2, 10), and the attenuation due to the channel is set equal to unity. The SINR 
is defined by SINR = 101og 10 [E\/(^ Y,k=2 E k + °"n))- In Figure the SINR is varied by changing the noise 
power a\ and the BEP is obtained for different SINR values in the cases of symbol-synchronous, chip-synchronous 
and asynchronous TH-IR systems with pulse-based polarity randomization and a synchronous TH-IR system without 
pulse-based polarity randomization 7 . For the asynchronous case, performance is simulated for different pulse shapes 
wi(t) and W2(t), given by d38l and d40l . respectively. From Figure |5j we see that the simulation results match 
closely with the theoretical results. Also note that for small SINR, all the systems perform quite similarly since 
the main source of error is the thermal noise in that case. As the SINR increases, i.e., as the MAI becomes the 
limiting factor, the systems start to perform differently. The asynchronous systems perform better than the chip- 
synchronous and symbol-synchronous cases since Jfr J C R 2 (e)de in d3Vj is about 0.2 for wi(t) and 2/3 for w 2 (t), 
which also explains the reason for the lowest bit error rate of the asynchronous system with UWB pulse wi(t). 
Also it is observed that for an IR-UWB system with pulse-based polarity randomization, the chip-synchronous and 

the symbol-synchronous systems perform the same as expected. Moreover, we observe that without pulse-based 

7 The results for the TH-IR system without pulse-based polarity randomization are provided to justify the discussion in Section III-D. The 
extensive comparison between TH-IR systems with and without polarity randomization is beyond the scope of this paper. 
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polarity randomization, the MAI is more effective, which results in larger BEP values. 

In order to compare the approximate analytical expressions and the simulation results for multipath channels, we 
consider the following channel coefficients for all users: a = [0.4653 0.5817 0.2327 - 0.4536 0.3490 0.2217 - 
0.1163 0.0233 — 0.0116 — 0.0023]. Then, the Rake combining fingers are (3 = ct for an ARake receiver, (3 = 
[0.4653 0.5817 -0.4536 0] for an SRake receiver with 3 fingers, and (3 = [0.4653 0.5817 0.2327 
0] for a PRake receiver with 3 fingers. The system parameters are chosen as N u = 10, N c = 5, 
N f = 15, E x = 0.5 and E k = 1 for k = 2, 10. Figure plots BEPs of different Rake receivers for synchronous 
and asynchronous systems with pulse-based polarity randomization. From the figure, we have the same conclusions 
as in the AWGN channel case about synchronous and asynchronous cases. Namely, chip-synchronous and symbol- 
synchronous systems perform the same and asynchronous systems with received pulses w\(t) and w^if) perform 
better. The asynchronous system with w\ (t) performs the best due to the properties of its correlation function. Note 
that the performance is poor when there is synchronism (chip or symbol level) among the users. However, the 
asynchronous system performs reasonably well even in this harsh multiuser environment. Hence, when computing 
the BEP of a system, the assumption of synchronism can result in over-estimating the BEP. Apart from those, it 
is also observed from the figure that the ARake receiver performs the best as expected. Also the SRake performs 
better than the PRake since the former collects more energy because the fourth path is stronger than the third path. 

For the next simulations, we model the channel coefficients as ai = sign(a/)|a/| for / = 1,...,L, where 
sign(cfy) is ±1 with equal probability and \a\\ is distributed lognormally as CM(fii,(J 2 ). Also the energy of the 
taps is exponentially decaying as E{|a/| 2 } = Qoe~ x ( l ~ l \ where A is the decay factor and X)^=i E{|az| 2 } = 1 
(so Qq = (1 — e _A )/(l — e~ Ai )). All the system parameters are the same as the previous case, except we have 
Ex = 1 in this case. For the channel parameters, we have L = 20, A = 0.25, a 2 = 1 and m can be calculated from 



m = 0.5 



ln( 1 i f^)-A(/-l)-2a 2 



, for I = 1,...,L. 

Figure plots the BEP versus E^/Nq for different Rake receivers in an asynchronous environment where wx(t) 
models the received UWB pulse. We consider ARake, SRake and PRake receivers for the TH-IR system with 
pulse-based polarity randomization and an ARake receiver for the one without pulse-based polarity randomization. 
The SRake and PRake receivers have 5 fingers each. As can be seen from the figure, the theoretical results are 
quite close to the simulation results. More accurate results can be obtained when the number of users is larger. It 
is also observed that the performance of the SRake receiver with 5 fingers is close to that of the ARake receiver in 
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Fig. 6. Bit error rate vs Eb/No for different cases, where N c — 5, Nf = 15, iV^ = 10, E\ — 0.5 and E = 1. The channel coefficients are 
[0.4653 0.5817 0.2327 - 0.4536 0.3490 0.2217 - 0.1163 0.0233 - 0.0116 - 0.0023]. 
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Eb/No (dB) 



Fig. 7. Bit error rate vs Eb/No for different receivers in an asynchronous environment, where N c = 5, Nf = 15, N u = 10, E\ = 1 and 
E = 1. The channel parameters are L = 20, A = 0.25, <r 2 = 1. The SRake and PRake have 5 fingers each. 



this setting. Moreover, the ARake receiver for the system without polarity randomization performs almost as worst 
as the PRake receiver for the UWB system with polarity randomization, which indicates the benefit of polarity 
randomization in reducing the effects of MAI. 

In Figure [8j we set E = 2 and keep all the other parameters the same as in the previous case. Here we consider 
a UWB system with polarity randomization and observe the performances of the SRake and the PRake receivers 
for different number of fingers M, using d29l . It is observed from the figure that the performance of the SRake 
receiver with 10 fingers is very close to that of the ARake receiver whereas the PRake receiver needs around 15 
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Fig. 8. Bit error rate vs Eb/No for different receivers in an asynchronous environment, where N c = 5, N f = 15, N u = 10, E\ = 2 and 
E = 1. The channel parameters are L = 20, A = 0.25, a 2 = 1. 

fingers for a similar performance. 

V. Conclusion 

In this paper, the performance of random TH-IR systems with pulse-based polarity randomization has been 
analyzed and approximate BEP expressions for various combining schemes of Rake receivers have been derived. 
Starting from the chip-synchronous case, we have analyzed the completely asynchronous case by modelling the 
latter by an equivalent chip-synchronous system with uniform timing jitter at interfering users. The effects of MAI 
and IFI have been investigated assuming the number of pulses per symbol is large, and approximate expressions 
for the BEP have been derived. Also for a large number of interferers with equal energy, an approximate BEP 
expression has been obtained. Simulation results agree with the theoretical analysis, justifying our approximate 
analysis for practical situations. 
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Appendix 

A. Asymptotic Distribution of n in @ 

The noise term n in © can be obtained from © and (Q as n = a n J s£| rtp (t)n(t)(it, where n{t) is a zero mean 
white Gaussian process with unit spectral density. Hence, n is a Gaussian random variable for a given template 
signal. Since the process has zero mean, n has zero mean for any template signal. The variance of n can be 
calculated as E{n 2 } = cr 2 J (s[]} mp (t)) 2 dt using the fact that n(t) is white. Using the expressions in Q and ©, 
we get 

(i+l)N f -l (i+l)N f ~l 
E{n 2 } = ^ £ J ff(t)dt + 2a 2 n £ df d^ J f 3 (t)f k (t)dt, (42) 

j=iNf j,k=iNf 

where f 3 {t) = £f =1 A w rx (t - jT f - cf ] T c - (I - 1)T C - n). 
It can be shown that / ff(t)dt = Ef=iA 2 for a11 

j since w rx (t) is assumed to be a unit energy pulse. Now 
consider f fj{t)f\{t)dt. By definition, fj(t)fk{t) is zero when there is no overlap between the pulses from the 
jth and the kth frames. Assume that L < N c . Then, fj(t)fk(t) = for \j — k\ > 1. In other words, there can be 
spill-over from one frame only to a neighboring frame. In this case, d42l becomes 




Note that fj(t) is a random variable at a given time instant t due to the presence of the random time-hopping 
sequence {cj }, and {fj{t)fj+i(t)} are identically distributed for j = iNf, + l)Nj — 2. Since {c^d^} has 
zero mean and forms an i.i.d. sequence for j = iNf, (i + l)Nf — 2, {d^d^^ J fj(t)fj + i(t)dt} forms a zero 
mean i.i.d. sequence. Hence, the second summation in (143 \ converge to zero as Nf — > oo, by the Strong Law of 
Large Numbers. 
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When the L < N c assumption is removed, we can still use the same approach to prove the result for finite values 
of L. In that case, we can write a more general version of ( 143 1 as 

L D (i+l)N f -l-k 

E{n 2 } = o*N f Y, $ + 2o* £ E / fj(t)f j+ k(t)dt, (44) 

Z=l fc=l i=iiVj. 

where fj(t)fk(t) = for |j — fe| > D. Since L is assumed to be finite, D is also finite. Hence, the second term in 
d44l) still converges to zero as Nf — ► oo. 

Thus for large Nf, E{n 2 } ff^iV/ X]^=i A 2 ' anc ^ so n * s approximately distributed as n ~ ^0 , o^iVy 

B. Proof of Lemma l3.il 

The aim is to approximate the distribution of a = xfj^Y^m^iN 1 ' 1 "m. where a m is given by (fTTTl . Note that 



N f £-tm=iN f 

a m denotes the interference to the mth frame coining from the other frames. Assuming that L < N c + 1, there can 
be interference to the mth frame only from the (m — l)th or (m + l)th frames. Hence, a m can be expressed as: 

a m = dW £ j ( (j - m)T f + - c ™) T • < 45 > 

j&{m— l,m+l} 

Note that , . . . , Qfi+iWy-i ^ identically distributed but not independent. However, they form a 1-dependent 
sequence [24] since d m and a n are independent whenever |m — n\ > 1. 

The expected value of a m is equal to zero due to the random polarity code. That is, E{a m } = 0. The variance 
of a m can be calculated from (145 \ as 



E{d 2 m } = E E | ((j - m)r/ + (c« - c«)T c )] 2 } 



(46) 



j'£{m— l,m+l} 

where the fact that the random polarity codes are zero mean and independent for different indices is employed. 

Since the TH sequence can take any value in {0, 1, . . . , N c — 1} with equal probability, the variance can be 
calculated as 

E ^ = ^2 T,3{[4^{JTc)] 2 + [4>$2HTe)] 2 } , (47) 
c j=l 
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which can be expressed as 



L-l 



'L-j 



(i) 



(48) 



a=i / m=i / 
using (fl"2l. © and ©. 

Now consider the correlation terms. Since L < N c , E{a m a n } = when \m — n\ > 1. Hence, we need to consider 
E{a m a m +i} only. Similar to the derivation of the variance, E{a m o m +i} can be obtained, from (145b . as follows: 



1 L-i /L-i \ /L-j 

E{d m a m+ i} = ^2 E ( E # a l& ) ( E ^ a ffi 



(49) 



j=l \l=l 



c- r- ■ , , , . [e7 ^{i+l)N f -l „ 

Since |am)m=i7v/ 1S a zero mean 1-dependent sequence, J jf- 2Jm =iNf a m converges to 



M(0,E 1 [E{(a m ) 2 } + 2E{a rn a m+1 }]) 



(50) 



as N f — ► oo [24]. Hence, O follows from (08j and (1491. 
C. Proof of Lemma \3.2\ 

In this section we derive the distribution of IFI for L > N c + 1. Consider the case where (D — l)iV c + 1 < 
L < DN C + 1, with D being a positive integer. Hence, {Sml^^N 1 f° rms a ^-dependent sequence in this case. 
Similar to Appendix |Bj we need to calculate the mean, the variance and the correlation terms for a m in (fTTT i. Due 
to the polarity codes, it is clear that E{a m } = 0. The variance can be expressed as follows, using (fTTTi and the fact 
that the polarity codes are zero mean and independent for different indices: 



E{a 2 m } = J2 E \ duv [U ~ m ) T f + ( c 



which can be calculated as 



j=-oo 



iVc-lJVo-1 oo 



= ^2 E E E [tfffi ((i - ^)^/ + (4 1} - c ™) T < 



i=0 Z=0 j=— oo 



(51) 



(52) 



using that fact that the TH sequence is uniformly distributed in {0, 1, . . . , N c — 1}. 
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Then, the variance term can be expressed as 



= W Ei { \*®UT e j\ a + \^(-jT c )} 2 } + kV0T c )l 2 + te(-iT c )l 2 ) , (53) 

iVc j=l 1 J c j=N c 1 } 



which can be obtained, using ill 2b . © and as follows: 



L~N a 



E{(a m ) 2 } = ^ E 



1 y 



7V2 



i=i 



vi=l 



+ (E«^+ 



vi=l 
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>i=l 



(54) 



Since djjV/ , • • • , 0(i+i)jv>-l f orm a ^-dependent sequence, we need to calculate E{a m a m+n } for n = 1, . . . , D. 
Then, the IFI term a in < fTUb can be approximated by 



D 



E{(ajAr / ) 2 } + 2 y~] E{ajjy f ajjy- f +n } 



n=l 



(55) 



as Nf — ► oo [24]. 

Using (fTTT i. (fT2l . © and ©, the correlation term in d55l) can be calculated, after some manipulation, as 



D 1 L—Nc / j \ / j \ 

EE{0iJV / a iW> +n} = E E^°i+i-i E a i ^-N-i 

n=l c j=l \i=l / \i=l / 

+ ^E^E^£j(E«^j- ( 56 ) 

Hence, (fl4l can be obtained by inserting (1541 and d5oT l into ( 1551 . 
D. Proof of Lemma \3.3\ 

In order to calculate the distribution of the MAI from user k, = \J~^Yl^r^iNf 1 a m\ we first calculate 
the mean and variance of a™ given by d!6t . where the delay of the user, r^, is an integer multiple of the chip 
interval: = AfcT c . 

Due to the polarity codes, the mean is equal to zero for any delay value r^; that is, Eja^^lA^} = 0. In order to 
calculate the variance, we make use of the facts that the polarity codes are independent for different user and frame 
indices, and that the TH sequence is uniformly distributed in {0, 1, . . . ,N C — 1}. Then, we obtain the following 



26 



expression: 



Na-1N C -1 OO 

'• ^0 1=0 j=-oo 



E{(aW) 2 |A fc } = ± g £ £ {^[^-^(i-H^c + A^rj} 5 

which is equal to 

E{(a«) 2 |A fc } = ± g [<£)(jT c 
Using (I12L © and ©, d58b can be expressed as 



(57) 



L-l 



j=-(L-l) 



(58) 



E{(a«) 2 |A fc } = 1 



L-l / j 



(59) 



user 



j=i \i=i / j=i \i=i / 

Moreover, we note that E{a m |Afcj = for m ^ n due to the polarity codes. 

Similar to the proofs in Appendices iBl and ICl {dm }^}^ f 1 forms a dependent sequence and the MAI from 
K oP^ = \J~^Y^rt^iN f 1 ~ 1 a m' converge to AT ^0 , £?fcE{(am') 2 }j since the correlation terms are zero. 
Hence, ( fT7t can be obtained from d59t . 

Note that the result is true for any value of A k since E{(am^) 2 | A^} in (1591 is independent of A k . Hence, the 
result is valid for both symbol and chip synchronous cases. 

E. Proof of Lemma 13.41 

The proof of Lemma l3~4l is an extension of that of Lemma 1331 Considering d22b . we have an additional offset 
e k , which causes a partial overlap between pulses from the template signal and those from the interfering signal. 

Due to the presence of random polarity codes, the mean of a™ in d22l is equal to zero. Using the fact that 
the polarity codes are zero mean and independent for different frame indices and that the TH codes are uniformly 
distributed in {0, 1, . . . ,N C — 1}, we can calculate the variance of a™ conditioned on A^ and e k as 



1 



7V C -1 AT-1 oo 



E{(aW) 2 |A fc , efc } = ^£ £ £ , 

c i=0 1=0 j=-oo 



[(i-l + U -m)N c + A k )T c + e k }} , 

- t=[) t=u j=-oo 

which can be shown to be equal to 

E{(a«) 2 |A fe ,e fe } = — £ (jT c + e k ) 



(60) 



i=-L 



(61) 
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Note that since the expression in ( 16 1 1 is independent of A k , E{(am^) 2 | A*., e/J = E{(a^) 2 |efc}. 
From (11 2t . © and ©, we can obtain an expression for 4>uv(jT c + e&) when j > as 

i-j'-i 



,0T C + e k ) = E a< i k) \Pi+iR( e k) + 0l+j+lR(T c ~ e k )} + af_^ L R{e k ), (62) 
i=i 



where R(x) = w rx (t — x)w rx {t)dt. Similarly, the expression for <p uv {—jT c + e k ) can be expressed as follows 
for j > 0: 

L-j 



4>uv{-jT c + e k ) = J2 A [a\%R(e h ) + af^_ x R(T c - e k )] + L ^ +1 af R(T C - e k ). (63) 
l=i 

Using ( 162t and (163b . E{(am^) 2 |efc} can be expressed from doTl as 

c j=0 \i=l J 
1 L-l / j \ 2 

+ AT E E a * (fc) + p i+L -iR{T c - e k )\ + af^ L R(e k ) . 



(64) 



3=0 \i=l 

,(*)„(*) 1 



Also, due to the polarity codes, the correlation terms are zero. That is, ^{a m a n } = for m / n. Then, from 
the central limit argument in [24], we see that in (fT3l . conditioned on e k , converge to the distribution given in 
Lemma 13.41 



F. Proof of Lemma 13.51 

Consider (N u — 1) interfering users, each with bit energy E. Then, the total MAI a = J2k=2 a ^ * s tne sum °f 
(N u — 1) i.i.d. random variables, where = \J~-§j Y^rv=iNj 1 a rn ■ Using the results in Appendix |E| namely, 
E{am'} = 0, E{om "n } = for m ^ n and (1541 . we obtain 

= E « (fc) ~ M ( , _EK M/ifc (6 fc )} , (65) 

! 1 V C / 



^ ~ k =2 



as iV u — ► oo, where E{cr 2 /Af k {e k )} can be obtained as in d28l from d64l using the fact that t k ~ £^[0, T^). 



